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ABSTRACT 

We report the first determination of the binary period and orbital ephemeris of the Be X-ray binary containing the pulsar 1 A 1 118-616 
(35 years after the discovery of the source). The orbital period is found to be Po,b = 24.0 ± 0.4 days. The source was observed by 
RXTE during its last large X-ray outburst in January 2009, which peaked at MJD 54845.4, by taking short observations every few 
days, covering an elapsed time comparable to the orbital period. Using the phase connection technique, pulse arrival time delays could 
be measured and an orbital solution determined. The data are consistent with a circular orbit, and the time of 90 degrees longitude 
was found to be T„/2 = MJD 54845.37(10), which is coincident with that of the peak X-ray flux. 
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1. Introduction 

The X-ray transient lA 1118-61 was discovered during an out- 
burst in 1974 by the Ariel-5 satellite in the 1.5-30keV range 
dEvles et al.lll975i) . The modula tion with a period of 6.75 min 
found in the data used by I ves et alj ( Il975h was i nitially inter- 
preted as the orbital period of two compact obiects. lFabian et aU 
then suggested that the observed period may be due 
to a slowly spinning accreting neutron star in a binary sys- 
tem. The optical c ounterpart was identified as t he Be-star He 
3-640/Wray 793 bv lChevalier & Ilovaiskvl (fT97l and classified 
as an 09.5IV-Ve star wi th strong Balmer em ission lines and 
an extended envelope by Janot-Pacheco et al. (1981). T he dis - 
tance was estimated to be 5 + 2 kpc ( Janot-P acheco et a l.iil981h. 
The cl assification and distance was confirmed bv ICoe & Pavnd 
(II985h using UV observations of the source. 

A second strong X-ray outburst occur red in January 1992, 
which was observed by CGi?0/BATSE (Coe etalJ[T994h . The 
measured peak flux was ~ 150 mCrab for the 20-100 keV energy 
, range, similar to that of the 1974 outburst. This outburst was fol- 
' lowed by enhanced X-ray activity for ~60 days, with two flare- 
like events about ~2 5 d and ~49 d after the peak of the outburst 
(see lCoe et al.1119941 Fig. 1). Pulsations with ~406.5s were de- 
tected up to lOOkeV and the pulse profile showed a single broad 
peak, close to a sinusoidal modulation with little change with 
photon energy. A spin-up with a rate of 0.016 s/day was observed 
during the decay of the outburst. Multi-wavelength observations 
revealed a strong correlation between t he Hn, equivalent width 
and the X-ray flux. This led ICoe et alj ( 11994 to conclude that 
expansion of the circumstellar disk of the optical companion is 
mainly responsible for the increased X-ray activity, including a 
large outburst if there is enough matter in the system. This con- 
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Fig. 1. The light curve of the outburst of lA 1118-616 as ob- 
served by RXTE /PC A (4^0 keV) in January 2009. The time 
resolution is 16 s. The continuous curve is the daily light curve 
as seen by RXTE/ASM (scaled to the PCA count rate and 
smoothed by taking the running mean of five consecutive days). 



elusion is supported b y pulsations from the source also being 
detected in quiescence (iRutledee et al.ll2007l) . 

The source remained quiescent until January 4, 2009 when 
a third outburs t was detected by Swift (Mangano et al. 200^ 
lManganol2009|). Pulsa t ions w ith aperiod of 407. 68+0. 02s were 
reported bv ^Man ganol (I2009D . The complete outburst was reg- 
ularly monitored by RX TE. INTEGRAL observed the source 
after the main outburst (iLevder et al.l '2009). Suzaku observed 
lA 1118-61 twice, once during the peak of the outburst and 
also ~20 days later when the flux returned to its previous level. 
lDoroshenkoetal.1 (l2010l) analyzed RXTE/PCA data of the most 
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Fig. 2. An example of a pulse profile (with 32 phase bins) of 
lA 1 118-616 as observed in -2.5 ks by RXTEfPCA (normalized 
to a mean flux of 1 .0), corresponding to the folding of 6 pulses 
with a pulse period of ~407.7 s. The center of the observing time 
is MJD 54844.375. The profile is repeated such that two phases 
are shown. This profile was used as a template in measuring rel- 
ative phase shifts of profiles from other integration intervals and 
the corresponding pulse arrival times. 



recent outburst with emphasis on spectral analysis, detecting a 
cyclotron resonant scattering feature (CRSF) at ~55 keV. 

In this letter, we report on an in-depth timing analysis of 
RXTE/PCA data over the entire outburst in January 2009, and 
the discovery of an orbital period of 24.0 d. This solution is sup- 
ported by the timing of the three large X-ray outbursts seen in 
history and by the analysis of pulse arrival times of the January 
1992 BATSE observations. 



2. Observations 

We used two data sets. First, and most importantly, we analyzed 
data obtained by RXTE/PCA during a 27 -day monitoring of 
lA 1118-616, which started on January 7, 2009, covering the 
complete outburst. The resulting light curve (4-40 ke V) is shown 
in Fig. [1] The source flux peaked on January 14 and decayed 
over ~15 days. The data structure is defined by the orbits of the 
RXTE satellite, with 29 RXTE pointings of a typical duration 
of around one hour. A total exposure of 87 ks was obtained. The 
PCA data were all taken in event mode ("Good Xenon"), pro- 
viding arrival times for individual events. The data were reduced 
using HEASOFT version 6.8. 

A second data set consists of pulse profiles generated from 
archival data of the observation of the January 1992 outburst by 
CGRO/BATSE, covering about 12 days. 



3. Timing analysis of pulse arrival times 

3.1. RXTE/PCA - 2009 

Some timing analysis of the R XTE data was performed ear- 
lier bv iDoroshenko et aP (l2010l) . who determined the pulse pe- 
riod and an initial value of its derivative but failed to identify 
the orbital modulation. Here, a more rigorous analysis is per- 
formed, which begins in a similar way by selecting 19 integra- 
tion intervals of the 29 pointings according to the time struc- 
ture of the individual satellite pointings (Fig.[T]i. These integra- 



Fig. 3. An example light curve with 16 s resolution of one con- 
tinuous observation by RXTE/PCA in the 4^0 keV interval, to- 
gether with the best-fit cosine function. 

Table 1. Orbital elements, Ppu^e and P^^ue of 1 A 1 118-616. The 
uncertainties are 1 cr (68%) for two parameters of interest. 



T,/2 [MJD (TDB)] 
Porb [d] 
a sin i [It-s] 
eccentricity e 
n [deg] 

Ppuisc [s] at Tjj.f' 

f pulse [SS~'] atTi-ef^ 



54845.37 ±0.10 
24.0 ' ± 0.4 3 
54.85 ± 1.4 
0.0 ' 
360.0 ' 

407.6546 ±0.0011 
(-1.8 ± 0.2) X 10-'' 



These parameters were kept fixed (see text), ^ Tef = 54841.259391 ' 
see text regarding this uncertainty 

tion intervals ranged from ~ 1500 s to ~ 15000 s, corresponding 
to between four and 38 pulses. Nineteen pulse profiles (with 32 
phase bins) were then constructed by epoch-folding barycenter- 
coiTected event times (since the Earth is always moving with 
respect to the source, the event arrival times are referenced to 
the center of the solar system). Fig. |2] shows one of these pulse 
profiles produced by accumulating six pulses (centered at MJD 
54844.375). Using this profile, the pulse arrival times (in MJD) 
for the other profiles were determined by template fitting. The 
uncertainties of these arrival times are of the order of 10 sec 
(~2.5 % of the pulse period). 

For the case of no binary modulation and non-zero first 
and second derivatives of the pulse period, the expected arrival 
times as a function of pulse number n are given by (see, e.g., 
iKellev et al.lfT980l: iNagasdl 19891: LStaubert et al.ll2009h 

f„ = fo + n + ^ «^ PsPs + \n^ pIp, + ■■■ (1) 
2 6 

With the given accuracy, it was possible to phase-connect the 
19 pulse arrival times and determine values for the pulse pe- 
riod and its first derivative by applying a quadratic fit in n. The 
sond derivative was not constrained and set to zero. In this fit, 
significant systematic residuals remained with a sine-like shape 
and an amplitude of several tens of seconds, indicative of pulse 
arrival-time delays due to orbital motion. Adding the additional 
term 



+ a sin / X cos[2;t (t - T„f2)/Poib] 



(2) 



yielded an acceptable fit as expected for a circular binary orbit, 
a sin / being the projected radius of the orbit in light-seconds. 
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Fig. 4. Delays of the pulse arrival times in lA 1118-616 for the 
outburst of January 2009 as observed by RXTE/PCA and best- 
fit sine curve for circular orbital motion with a period of 24.0 d. 
The dots around zero are the residuals to the best-fit solution. 



and T;r/2 (= T90) in MJD being the time at which the mean 
orbital longitude of the neutron star is 90°, corresponding to 
the maximum delay in pulse arrival time. We found that a sin / 
~ 55 light-s, Porb ~ 24 d, Tgo MJD ~ 54845, Ppuise ~ 407.7 s, and 
Ps ~ -1.9x lO^^ss '. 

We then realized that by using the original light curve, a sam- 
ple of which is shown in Fig. [3] we were able to determine the 
pulse arrival times more accurately, when fitted piecewise with a 
cosine function (keeping the pulse period fixed at the previously 
found best-fit values, and leaving the mean flux, the amplitude, 
and the zero phase as free parameters). The same 19 integra- 
tion intervals as before were used. The phase connection analy- 
sis of the pulse arrival times determined in this way then leads to 
the orbital elements and pulse period information summarized 
in Table [1] The pulse delay times and the residuals against the 
expected delays for the best-fit circular orbit is shown in Fig. |4] 
The;\f^ of this best fit was 21.0 for 15 dof (degrees of freedom). 

In finding this solution, we assumed a vanishing eccentric- 
ity and a fixed value of 24.0 d for the orbital period. When the 
orbital period was left as a free parameter, we found that Porb - 
24.0 + 0.4 d. This uncertainty is expected to be large since our 
observations cover only 27 days, little more than one orbital cy- 
cle. 

Finally, we attempted a proper orbital solution with the ec- 
centricity £ and the longitude of periastron passage Q included as 
free parameters by solving Kepler's equation. The data do allow 
one to define a pair of parameters with reasonably constrained 
uncertainties of e = 0.10 + 0.02 and = (310 + 30) °. However, 
we consider the evidence of a non-zero eccentricity as marginal. 
By introducing these two additional free parameters, the was 
reduced from 2 1 . ( 1 5 dof) to 1 6 .4 for ( 1 3 dof) . An F- test yielded 
a rather large probability of ~22 % for the improvement to have 
occured by chance. Furthermore, no systematic uncertainties of 
any kind were considered. 

We note that the values for the pulse period and its derivative 
given in Table 1 differ from those stated by Doroshenko et al. 
3201 0) because of the missing orbital corrections in the earlier 
analysis. 



Fig. 5. Delays of the pulse arrival times in lA 1118-616 for the 
outburst in January 1992 as observed by CGRO/B ATSE. 

3.2. CGRO - 1992 

The second data set is from observations of the January 
1992 o utburst by CGRO/BATSE. This outburst is described 
well in ICoe et al.l (1994). Pulse profiles were generated using 
phase/energy channel matrices for all eight BATSE detectors 
which are stored at the HEASARC archive. Q For each detec- 
tor, events in the 20-40 keV energy range (according to the en- 
ergy calibration provided) were selected and sorted into com- 
mon pulse profiles with 128 phase bins. This was done for 12 
different integration intervals, covering about 12 days of obser- 
vation (MJD 48621^8633). Pulse arrival times in MJD were 
determined by template fitting. Fig. |5] shows the delays of pulse 
arrival times as found from a pulse phase connection analysis. 
Assuming an orbital period of 24.0 d, the best-fit solution for a 
circular orbit leads to the parameters of Ppuise - 406.53 + 0.02 s, 
/"pulse = (-3.1 +0.9)x 10-^ss-i, and T„/2 = MJD 48633.5 +2.5 d. 
When the difference between the Tjr/2 value from the RXTE ob- 
servation (see Table 1) and the one from the BATSE observation 
is divided by the orbital period of 24.0 d, we find a separation 
of 258.83 orbital cycles. If, in turn, we divide the separation by 
259 cycles, we determine a cycle length of 23.98 d. The corre- 
sponding uncertainty in this value is 0.01 d. However, we cannot 
be certain wether the separation is really 259 cycles, since the 
previously found uncertainty of 0.4 d in the orbital period would 
permit any cycle number between 255 and 263. 

4. Other support for the 24.0 d orbital period 

4.1. Timing of large X-ray outbursts 

Three large X-ray outbursts of lA 1 118-616 have been observed 
so far. The first one occurred in December 19 74, leading to the 
original discovery of the source by ARIEL-5 (' Evles et al.l[T975t 
Ives et al. 1975). At this time, a modulation with a period of 
6.75 min was also discovered. The s econd burst wa s observed 
by CGi?0/B ATSE in Januai-y 1992 (iCoe et al.lll994i) . with en- 
hanced X-ray activity for ~80 days. The third, most recent out- 
burst occurred in January 2009 and was observed by several 
X-ray sa teflites: Swift (M angano et al. 2009; Mangano 200^, 
RXTE (iDoroshenko et all 120101) . INTEGRAL dLevder et d] 

' ftp: //heasarc. gsfc.nasa.gov/compton/data/batse/pulsar/ 
ground folded/All 18-6 
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Fig. 6. A section of the X-ray light curve of lA 1 118-616 as ob- 
served by the RXTE/ ASM (smoothed using a running mean of 
five consecutive days). The large burst of January 2009, reach- 
ing a peak flux of 7.6cts/s at MJD 54845.4 is included. The 
short vertical lines are at phase 0.0 with respect to the ephemeris 
with Porb - 24.012 d and the arrows point to those small flares 
with peak fluxes >0.5 cts/s (dotted horizontal line), which hap- 
pen close to phase 0.0 (several others appear to happen close to 
phase 0.5). 



l2009 l). and Suzaku. Also in this case, the source continued to 
exhibit high activity for about 70 days after the large outburst. 
From the cited publications, we determined the time of the peak 
fluxes in these outbursts to be MJD 42407.0, MJD 48626.0, 
and MJD 54845.4, respectively, with an estimated uncertainty 
of +1 day in all cases. Taking the above value of 24.0 d for the 
orbital period, the outbursts in December 1974 and the one in 
January 2009 occurred 259 orbits before and after the one in 
January 1992. But, again, because of the uncertainty in Porb the 
separations could range from 255 to 263 orbital cycles, corre- 
sponding to a series of period values separated by about 0. 1 d. In 
Section 4.3, however, we present evidence that 259 may be the 
correct number Assuming that this is indeed so, a linear fit to 
the three outburst times leads to a period of (24.012 + 0.003) d 
(the small uncertainty being due to the long baseline in time). 

4.2. Timing of medium-size X-ray flares 

As already mentioned, the source remains particularly active af- 
ter the second as well as after the third l arge outburst for several 
tens of days. Fig. 1 in ICoe et"al ] (11994 shows this for the sec- 
ond large burst (observed by CGRO/BATSE in January 1992): 
there is highly structured X-ray emission with several peaks, the 
largest of which are around ~26 d and ~49 d after the peak of 
the large burst (with peak fluxes ~0.5 and ~0.4 of the big burst, 
respectively). In Fig. |6l we show the light curve in the vicin- 
ity of the third large burst (of January 2009) as observed by 
J?XTE/ASM (in a smoothed version of the daily light curve with 
a running mean of five consecutive days). Three peaks can also 
be identified within ~70 d of the the large burst, none of which 
corresponds exactly to phase zero, but the second and the third 
one are close (the mean of the three separations - starting with 
the large burst - is ~23 d). 
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Fig. 7. Frequency histogram of ASM small flares with peak flux 
>0.5 cts/s (from smoothed daily light curves) as a function of the 
24.012 d phase. 



4.3. Timing of small X-ray flares 

In the section of the ASM light curve shown in Fig.|6l three small 
peaks before the big burst (highlighted by small arrows) happen 
quite close to phase zero. There are also flares at other phases, 
but counting the number of flares with peak fluxes higher than 
0.5 cts/s over the complete ASM light curve (MJD 50087 - MJD 
55315), we find 39 out of 86 flares (45 %) between phase -0.15 
and -hO. 15. The rest are almost uniformly distributed with a slight 
enhancement around phase 0.5. The analysis of excursions with 
fluxes > 1 .0 cts/s confirms this result: 15 out of 32 events between 
phase -0.15 and -h0.15. Fig. Q shows the frequency distribution 
for the case of the >0.5 cts/s peaks. When the same analysis is 
repeated with orbital periods corresponding to separations be- 
tween the large outbursts other than 259 orbital cycles (that is 
between 255 and 263), the rate of coincidences of small flares 
with phase zero is significantly less. We take this as an indica- 
tion, not proof, that 259 is probably the right number. 

When a simple epoch folding of the long-term light curve 
of 1 A 1 118-616 as observed by RXTE/ASM is done, no signif- 
icant peak is found at ~24d. This is not inconsistent with the 
above considerations of medium size and small X-ray flares: 
there are only very few (not well aligned) medium size flares, 
occurring after the big outbursts, and the small flares (>0.5 cts/s) 
have such low fluxes that both types of flares are completely 
buried in the "noise" represented by the daily flux measurements 
<0.5 cts /s. When a dynam ical power density spectrum (PDS) 
(see e.g.l Wilms etal]|2001h is generated (with an individual data 
set length of 2000 d and a step size of 10 d), a weak signal in 
the form of a broad peak between 22 d and 25 d is found. The 
PDS is capable of detecting quasi-periodic signals allowing for 
frequency variations with time. 

5. Inclination and mass function 

Using the orbital elements determined above and an estimate of 
the mass of the optical companion, we can constrain the inclina- 
tion of the binary orbit o f lA 1118/He 3-640 . The o ptical com- 
panion was classified by Janot-Pacheco et alj (Il98lh as an 09.5 
Ill-Ve sta r with strong Balmer emission lines and an extended 
envelope. iMotch et al.l (11988) assumed a mass of 18 Mq for the 
optical companion. When the calibration of O-star parameters of 
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iMartins et alJ (l2005h is used, the value of 1 8 M© is confirmed (in- 
terpo lating the values for 09.5 of Tables 4 and 5 of' Martins et al.l 
l2005h . Adopting this value and assuming 1.4 Mq for the mass of 
the neutron star, Kepler's third law leads to a physical radius of 
the (circular) orbit of a = 219.1 light-s. The observed a sin i = 
54.85 light-s then leads to sin i = 0.25 and to the inclination of 
i - 14.5°. This value is consistent with no eclipses being ob- 
served. The mass function is f(M) = (m2 sin i)-' / (mi-i-m2)^ - 
1.14 X 10-1 

6. Discussion 

The detection of the period of 24 days for the binary orbit of 
1 A 1 118-616 rests mainly on the observation and analysis of or- 
bital delays of the arrival times of the X-ray pulses (with a pulse 
period of ~407 s). The main data set is from observations by 
RXTE/PCA, which sampled the large X-ray outburst of January 
2009. Fortunately, the duration of the outburst as well as the 
length of the observations were long enough to cover slightly 
more than one complete orbit. It was also fortunate that the sam- 
pling pattern was dense enough to allow pulse phase connection 
between the sampling intervals. However, with data for only one 
orbit the determination of the orbital period using these data has 
a rather large uncertainty of about +0.4 d. We asume this value 
of 0.4 d as the final uncertainty in the orbital period, despite the 
evidence from the small flares that 24.012 d, corresponding to 
a separation of 259 orbital cycles between the three large out- 
bursts, may be the correct value. 

For the January 2009 outburst, we have found that the peak 
of the X-ray flux coincides with T„/2, i.e., an orbital longitude 
of 90 °, while the formal value of O was determined to be 310 °. 
This supports our reservation of taking the formally achieved 
e / Q combination seriously: an F-test does give a rather high 
probability (~22 %) that the improvement in (when these 
two free parameters are introduced) is just by chance. On the 
other hand, the finding that a majority of the small flares and 
one of the large bursts do occur at or around phase zero of our 
24.012 d ephemeris indicates that the eccentricity may be some- 
what larger than zero, but most likely <0.16. 

For completeness, we mention here that in optical observa- 
tions of He 3-640/Wray 79 3 large values of l i„ equivalent width 
are occasionally observed. ICoe et aT l ([T99l had found a value 
in excess of 100 Angstrom on MJD 48706 (80 d after the peak 
of the January 1992 large X-ray outburst) and they associated 
the two phenomena with each other, which is clearly justified 
by examples of th ese associations in other B e X-ray binarie s 
dGrundstrom et al. !2007: Kiziloglu et al. 2009; Reig et alJ20"Toh . 

With respect to the timing of the three observed large X-ray 
bursts, w e emphasize t hat th e two separations are the same at 
17.04 vrs. lMotch et al.1 (Il988l) already concluded that the enve- 
lope of He 3-640 is probably not in a stationary state but under- 
goes expansion and contraction phases on a time scale of several 
years." We suggest that the 17yrs between the large X-ray out- 
bursts is a characteristic period for the oscillation of the envelope 
of He 3-640. 

Finally, we look at the position of lA 1118-616 in the 
Corbet-diagram, wh ich relates the o rbital p eriod to the spin pe- 
riod (see e.g.. iReig et al]l2004t iRodriguez et al.ll2009l) . For Be 
binaries, there is considerable s catter around a mean correlation 
trend, which was quantified by ICorbell (Il986h by the following 
formula: = 10 days x (1 - e)-2/3 (p^^.^ji /)i/2 According 
to this formul a, an orbital period of ~200d or more had been 
expected (e.g.. lMotch et al.lll988h . a prediction that in addition 
to the generally low level of the X-ray flux and the rareness and 



shortness of the larger outbursts may have helped the 24 d orbital 
period escape detection for 35 yrs after the source's discovery in 
1975. lA 1118-616 is indeed at the edge of the Be star distri- 
bution towards the region where wind-fed SGXRBs (super giant 
X-ray binaries) accumulate, and the sec ond most off'set Be star 
after SAX J2103.5-H4545. For the latter. iReig et"al] (l200l had 
discussed a physical reason for the slow spin at a rather short 
orbital period, namely that long episodes of quiescence in Be 
X-ray systems can cause a spin down to an equilibrium period 
that is expected for wind-driven accretion. This may also apply 
to lA 1118-616 and is consistent with the observed moderate 
spin-up rate (~1.1 Hz s"') during its outbursts in January 2009. 
A discussion of variou s alternative w ays of reaching very long 
spin periods is given in lFarrell et al.l (j2008) in the context of in- 
vestigating the extremely slow (~2.7 h) pulsar 2S 01 14-1-650. 
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